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(as  used  in  [14])  but  analogous  to  lrreduclblllty 
In  Markov  chain  theory,  and  which  guarantees  that 
the  invariant  measure  has  a certain  Isotropic 
property  for  ||x||  large. 


In  this  paper  we  consider  the  behavior  for 
large  t of  the  solutions  of  Ito  equations  of  the 
form 

dx  - Axdt+EB.xdw  +£g  dw 


2.  Motivation 

We  recall  here  a few  examples  from  [14]  and 
some  comments.  For  the  scalar  Ito  equation 

dx  * axdt+bxdWj+cch^ 
there  is  an  invariant  measure  given  by 


In  particular  we  investigate  properties  of  the  In' 
variant  measures  of  such  systems  and  describe  a 
certain  Isotropic  characteristic  of  the  Invariant 
measure. 
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provided  c ^ 0 and  a < 0.  We  see  In  this  case  that 
only  a finite  number  of  moments  exist  for  the  limit- 
ing distribution,  regardless  of  how  many  exist 
Initially.  By  Juxtaposing  two  or  more  such  equations 
with  different  values  of  6 we  can  construct  In- 
variant measures  for  bilinear  systems  which  have 
the  form 


Introduction 

The  steady  stream  of  papers  on  bilinear 
systems  In  recent  years  is  due  in  part  to  the  very 
general  nature  of  this  class  of  systems  and  in  part 
to  the  rather  satisfying  degree  to  which  they  can 
be  analyzed.  For  example  there  are  now  two 
different  kinds  of  approximation  theorems  [1,2] 
and  [3,4]  which  Illustrate  the  approximation  of 
quite  general  nonlinear  systems  with  bilinear  ones 
and  there  is  a general  Isomorphism  theorem  for  bi- 
linear systems  [5,6,7]  which  uses  controllability 
and  observability  to  prove  a uniqueness  theorem 
for  bilinear  realizations.  The  Volterra  series  is 
a perfectly  satisfactory  description  of  the  input- 
output  map  defined  by  bilinear  systems  [8]  and  the 
realization  problem  for  bilinear  systems  is  now 
reasonably  well  understood  [9,10,11].  In  our 
earlier  papers  [12,13,14]  we  have  undertaken  a 
study  of  bilinear  stochastic  models.  In  the 
stochastic  case  the  problems  seem  to  be  more 
difficult  and,  with  the  exception  of  those  results 
which  have  been  obtained  by  analysis  of  the  moments, 
rather  little  Is  known. 

In  this  paper  we  continue  our  earlier  study 
of  the  limiting  probability  distribution  associated 
with  Ito  equations  of  the  form 
m r 

dx  - Axdt  + l B.xdw.  + l g.dw.  (*) 

>-  i-1  x 1 i-m4-l  1 1 

~~  Our  main  results  concern  rather  precise  estimates 
Cj  on  the  asyq>totic  decay  (for  large  ||x||)  of  the 
CJ>  invariant  measures.  We  introduce  a new  nixing 
condition  which  Is  stronger  than  controllability 


( ox  x^+1 ) Bl  (o2  x*+l ) . . . (anxj>l ) 6x1 
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These  systems  are  "fully  elliptic"  In  that  there 
are  as  many  (in  fact  twice  as  many)  independent 
Wiener  processes  as  dim  x. 

On  the  other  hand,  we  show  in  [14]  that  the 
bilinear  model  (*)  can  have  a gausslan  invariant 
measure.  In  particular,  for  the  Ito  equation 

dx  • (S-bb')xdt-Mlxdw.+bdw, 


SI  and  SI  both  skew  symmetric  the  gausslan 


Is  Invariant.  This  construction  Is  valid  for  x 
of  any  dimension,  assuming  controllability,  and 
hence  may  describe  a highly  degenerate  elliptic 
situation. 

This  then  leads  to  questions  about  the  "typical" 
asymptotic  behavior  of  bilinear  stochastic  equations. 
In  particular,  we  are  Interested  in  the  following 
question.  Is  it  to  be  expected  that  the  existence 
of  moments  is  isotropic,  l.e.  independent  of  the 
direction  in  x-space,  as  it  is  for  the  second 
example,  or  is  It  to  be  expected  that  a different 
niaber  of  moments  will  exist  depending  on  the 
direction  as  in  the  first  exanpleT  Our  results 
here  show  that  the  first  alternative  is  the  rule 
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and  the  second  Is  In  some  sense  exceptional. 

3.  Moment  Observability 

The  main  tools  available  so  far  to  study  bi- 
linear stochastic  equations  are  the  associated 
moment  equations  and  certain  Ideas  from  controll- 
ability and  hypoelliptlclty . We  recall  from  [14] 
that  If  we  Introduce 

_P  1 
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where  the  meaning  of  (p ] Is  to  be  derived  from 
expressing  the  derivative  of  xIpJ  along  solutions 
of  x ■ Fx  as 


. d [p  1 
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(See  [11-14]  for  more  details.)  The  matrix  G Is 
constructed  in  [14].  p 

We  recall  that  for  A a square  matrix  and  c 
a vector,  the  pair  (A,c)  is  called  observable  if 
c is  a cyclic  vector  for  A’  . This  Is  necessary 
and  sufficient  for  the  initial  state  of  the  free 
dynamical  system 

x - Ax;  y - <c,y> 

to  be  determinable  from  the  observation  of  y and 
It  Is  sufficient  to  insure  that  y(-)  Is  boixided 
on  [0,«)  If  and  only  if  eAtx(0)  is  bounded  on 
I0,«). 

In  two  dlmenslona  it  happens  that  for  some 
matrices  A,  any  c i 0 gives  an  observable  pair. 

More  precisely,  in  this  case  if  A has  complex 
eigenvalues  and  c is  nonzero  then  (A,c)  is  an 
observable  pair  (we  are  assisting  A and  c are  de- 
fined over  the  real  field).  However,  for  dimen- 
sions greater  than  2 there  is  for  each  A,  vectors 
c(A)  such  that  c(A)  « not  zero  and  yet  (A,c(A))  is 
not  an  observable  pair.  We  make  these  coassnts  by 
way  of  contrast  to  what  follows. 

Consider  the  Ito  equation  (e)  and  the 
associated  moment  equation  (t).  We  focue  atten- 
tion on  the  homogeneous  part  of  the  moment  aquation. 
The  Ito  equation  ( * ) will  be  called  p- moment 
observable  If  for  every  c d 0 the  pair 

((A_  T M[>(p]  +Tr(,i(P])2'  'IP,> 


la  observable.  What  this  definition  means  la  this. 
Mote  that  (sea  [12]) 

so  that 

#<c,x>p  - <c<”. 


thus  we  are  isolating  the  condition  under  which  we 
can  determine  any  pth  moment  from  the  time  history 
of  any  other  pth  moment.  Only  observational  vectors 
of  the  form  clp),  a rather  thin  subset  of  all 

( p ) vectors,  have  meaning  in  this  context. 

Recall  that  a set  of  square  matrices  ^BjJ  Is 
said  to  act  lrreduclbly  if  there  is  no  cnange  of 
basis  such  that  they  all  take  the  form 


PBjP-1  - 


"11 


B12 

B22 


It  should  be  clear  that  if  {B. } is  reducible  then 


(B 


Up] 


} is  alsoj  as  is  the  whole  associative 


algebra  generated  by  the  Bf's  or  the  B.  , ,'s. 

Thus  we  see  that  the  following  lemna  is  Valid. 

Lenma  1:  A necessary  condition  for  p-moment 
observability  of  (*)  is  that  (A,B^ .B^ . . .B^} 
should  act  lrreduclbly. 

In  order  to  develop  additional  intuition  about 
the  question  of  p-moment  observability  we  point 
out  that  for  (A,b,c)  of  appropriate  dimensions 


(<c,eAtb>)P  - <c[P',eA(Pjtb(PJ> 


but  that  if 

<c,eAtb> 
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,b^,c^)  is  controllable  and  observable 


Thus  if  the  (n+£  X)  in  number  p-fold  sums  of  the 
eigenvalues  are  distinct  then  the  triple 

(Alp] 

by  well  known  linear  system  theoretic  arguments. 
This  leads  to  the  following  conclusion. 

If  A has  distinct  eigenvalues  say 

and  if  the  p-fold  sins  of  these 


X1,X2 Xn 

eigenvalues  are  distinct  then  b is  a cyclic  vector 
for  A if  and  only  if  b IP J is  a cyclic  vector  for 
A,_t.  Likewise  c is  a cyclic  vector  for  A ' 


'(pr 


if  and  only  if  c^  la  a cyclic  vector  for  Aj^j. 

We  now  address  the  question  of  existence  of 
p-moment  observable  systems. 

Thsorem  1 : The  stochastic  equation 

dx  * Axdt+BxdWj+gdWj  (**) 

is  generlcally  p-moment  observable. 

Remark i What  we  mean  by  this  is  that  there  is  an 
open  and  dense  aidteet  of  7fnan  » gjn*n  auch  that 
for  A and  I in  this  subset  we  have  p-moment 
observability. 

Proof;  The  sat  of  A'a  which  have  no  repasted 
eigenvalues  is  an  opan  and  danse  set.  If  the 
eigenvalues  of  A are  isi  rapes  ted  we  can  put  it  in 


diagonal  form. 


In  this  case  is  also  diagonal. 
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The  Invariant  subspaces  of  A In  this  case  are  Just 
the  subspaces  characterized  by  vectors  with  a 
certain  number  of  zero  cong>onents.  There  are,  of 
course  2n  such  subspaces  corresponding  to  all 
possible  ways  of  fixing  zero  elements.  Now  con- 
sider the  moment  equations 

Tt'*[p]  ‘ «A- 1®2) iPi+  i<BlP])2,,r,tlpl+Gp  *lp'21 

Ue  let  A have  unrepeated  eigenvalues  and  be  In 
diagonal  form.  We  are  Interested  In  the  set  of 
vectors 


l(A‘  Ib2)[p1+I(b1p1)2>,,Cc1p1: 


Now  the  determinant  of  matrix  formed  by  these 
columns  can,  for  each  value  of  A and  c,  be  regarded 
as  a polynomial  in  the  entries  of  B.  We  want  to 
show  that  this  polynomial  does  not  vanish  identi- 
cally ixiless  c Is  zero.  From  our  previous  remarks 
we  see  that  the  zeroth  degree  term  In  the  polynom- 
ial vanishes  only  if  c has  one  or  more  zero  entries. 


By  reordering  the  variables  we  can  arrange 
matters  so  that  a given  invariant  subspace  of  A 
consists  of  all  vectors  of  the  form 


S(A)  - 


S(A[p]>- 


}r- components 
} (n-  r)- componen  ts 

}(n+P  ^)_(n  T"P  *)-cong>onent8 
P P 

)(n  r~P”1)-coiqponents 


Now  for  this  subspace  to  remain  invariant  when  we 
12  1 2 

add  jfl  | ^(Bj  j)  ft  is  necessary  that 


But  for  generic  B no  such  zero  block  exists  and  so 
generlcally  we  have  p-aoment  observability. 

4.  Aaynitotlc  Properties  of  Invariant  Measures 

We  now  apply  the  previous  results  to  the 
study  of  the  possible  forms  of  Invariant  measures 
achlewd  as  the  limiting  distribution  of  stochastic 
differential  equations  of  the  form  here.  The  main 
thrust  of  the  following  theorem  together  with 
Theorem  1,  is  to  show  that  typically  when  limiting 
distributions  exist  they  aim  isotropic  In  the  sense 
that  their  rata  of  decay  along  any  line  In  x-space 
is  roughly  the  sans . 

Theorem  2;  If  the  moments  of  (*)  up  to  and  in- 
cluding order  p-1  are  asymptotically  stable  and  if 
the  monant  aquation  of  order  p la  instable  then  for 
no  c does 

111  #<<c,x(t)>)p 

t-*» 

exist,  provided  (*)  la  p- moment  obaervable. 

Proof:  This  la  an  easy  consequence  of  the  earlier 
remarks  on  observability  and  the  da  fin  it  loo  of 
p- no  sent  observability. 


One  reason  this  result  seems  Interesting  Is 
that  it  points  out  more  precisely  what  forms  the 
invariant  measures  for  bilinear  stochastic  systems 
can  take.  In  particular  It  can  be  used  as  a guide 
In  attempting  to  find  series  expansions  for  the 
Invariant  measures. 
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